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A cluster growth model is proposed to study chiral symmetry breaking in stirred crystal- 
lization from an achiral element. Achiral monomers are assumed to coagulate to form chiral 
clusters from dimers to hcxamcrs. Due to the stirring, the hexamers break into dimers. The 
coagulation of two dimers into a tetramer also occurs. The fixed point analysis of coupled 
rate equations of cluster densities and their numerical integrations show that the chiral sym- 
metry becomes broken if dimers are critical and able to dissociate back to achiral monomers. 
The chirality selection is understood by showing that, in a quasi-steady approximation, the 
rate equations reduce to those of dimers with nonlinear autocatalysis and decomposition. 

KEYWORDS: homochirality, crystal growth, enantiomer, rate equation, dynamical chiral sym- 
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It has been known for a long time that organic molecules in life are homochiral with a 
completely broken chiral symmetry, and their origin has intrigued many scientists. 1 Frank 
showed theoretically that an autocatalytic reaction with cross inhibition between two ena- 
tiomers leads to the selection of one type of enantiomer in an open system. 2 We have recently 
proven that a system with nonlinear autocatalysis and back reaction as is described by the 
following rate equations 

r = {ko + k\r + k2r 2 )a — Xr, 

s = (ko + k±s + k2S 2 )a — Xs, (1) 

can lead to complete homochirality in a closed system. 3-5 Here, a denotes the concentration 
of an achiral substrate, r and s are those of chiral enantiomers, and the conservation relation 
a + r + s =const. is assumed to hold. The model is motivated by the chemical reaction 
system discovered by Soai et al.®' 7 which is the first example of chemical reaction to show 
the amplification of enantiomeric excess in a closed system. 

Chiral symmetry breaking, however, is found not only in biological and chemical systems, 
but also in other systems such as in crystallization. 8 Even achiral inorganic molecules can 
gather together to form chiral crystals. Examples are a quartz of SiC-2 crystal, or a crystal 
of sodium chlorate (NaC10,3). When the sodium chlorate crystal grows from a solution in a 
cell, a racemic mixture of D and L crystals emerges. On the other hand, when the crystal is 
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grown from a supersaturated solution under continuous stirring, the growth cell is occupied 
by one type of enantiomer. 9 The chirality selection is found to be brought about by frag- 
menting an initially nucleated crystal: the scattered clusters of small crystal act as centers of 
secondary nucleation of the same chirality as that of the initial crystal. This senario is called 
a " single mother" crystallization. Recently, however, it has been found that even if the initial 
solution contains crystalline seeds of both enantiomers, the continuous stirring leads to the 
state of complete chiral purity with one type of enantiomer. The selected enantiomer is the 
one which was in a slight excess initially. 10 There is a theoretical study associated with this 
phenomenon, 11 but we propose another approach in order to examine the roles of various 
elementary processes. We also study the relationship between the present model and eq. (1) 
reported in our previous studies. 3,5 

Our model of crystallization is described in terms of clusters of various sizes. An individual 
molecule A is achiral, but a few of them coagulate into a chiral embryo. For simplicity, we 
assume here that dimers already show the chirality. We call the right-handed dimer R2, and 
the left-handed one S2- The first growth process is then described by 

A + A^R 2 , A + A%S 2 . (2) 

A A 

These chiral embryos grow further by incoorporating elements A one by one as, 

k- k- 

Ri + A-^Ri+x, Si + A-^Si + x, (3) 

with i > 2. Here, we assume for simplicity that the critical size of the crystal is i = 2. Only 
the dimers can dissociate with a rate A, which is assumed larger than the growth rate kx 
from monomers into a dimer. Clusters larger than the dimer are stable and only grow. If 
these are the whole processes, the crystalline enantiomeric excess does not increase. In the 
experiment 10 the solution is continuously stirred by a stirrer and balls, 10 which crash part of 
the grown crystals and split off small clusters. 

First we consider the case that the smallest nontrivial cluster, tetramer, splits as 

Rx^R 2 + R 2 , Si^S 2 + S 2 (4) 

with a rate v. Coupled rate equations for concentrations n, Sj of dimers (i = 2), trimers 
(i = 3) and tetramers (i = 4) are written as 

r 2 = k\c? — \r 2 — k 2 r 2 a + 2ur^ 

h = k 2 r 2 a - k 3 r 3 a 

^4 = k 3 r 3 a — vt/x (5) 

and the corresponding ones for S's. Since the total amount of the molecules is conserved in a 
closed reaction vessel, the concentration of elementary molecule A is given as a = c$ — 2{r 2 + 
s 2) — 3(r3 + S3) — 4(r4 + S4) with a constant cq. Numerical integration indicates readily that the 
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chiral symmetry will not be broken in the present case. This result can be easily undestood by 
applying the quasi-steady state approximation. Because of large reaction rates k{ (i > 2) and 
v, clusters larger than the dimers are assumed to adjust their concentrations rapidly to their 
steady-state values as r% = (k 2 /k 3 )r 2 , r 4 = {k 2 /v)r 2 a, as long as a 7^ 0. Then the evolution 
of the dimer concentrations, r 2 and s 2 , is simply described by the reduced rate equations 

f 2 = kia 2 - \r 2 + k 2 ar 2 , 

s 2 = kia 2 - \s 2 + k 2 as 2 . (6) 

The equation contains only a linear autocatalytic process and a back reaction, and looks similar 
to eq. (1) with k 2 = 0. 3 A simple calculation shows that there is only a racemic fixed point as 
r 2 = $2 = k\a? /(A — k 2 a). Therefore, the system approaches a racemic state asymptotically, 
in agreement with numerical integration results. The chiral symmetry is conserved. The key 
point is that the secondary nucleation induces the autocatalytic effect in the dimer evolution, 
but only up to the level of linear autocatalysis. Note that the fragmentation effect v does not 
appear explicitly in the reduced rate equation (6). 

Since the chirality is not selected in the cluster growth up to the tetramer, we go one step 
further, up to the hexamer; Rj, Sj with 2 < i < 6. Instead of eq. (4), we assume the breaking 
of a hexamer into three dimers as 

Rq A 3R 2 , Sq —* 3S 2 , (7) 

and furthermore the incorporation of the coagulation of dimers into a tetramer 

R 2 + i?2 — ► i?4, 5*2 + 1S2 — ► S4 (8) 

Then the rate equations are written as 

r 2 = k\a 2 — \r 2 — k 2 r 2 a — 2fir 2 + 3vr$ 

h = k 2 r 2 a - k 3 r 3 a 

^4 = A;3T3a — k^r^a + p,r\ 

r5 = k^r^a — k^r^a 

h = k 5 r 5 a - ur G (9) 

and the corresponding ones for S's, under the conservation condition a = Co — Y^i=2 ^( r * + 

Since these coupled equations involve many degrees of freedom, detailed analysis can 
be performed only by numerical integration. However, the essential feature may be under- 
stood in a quasi-steady approximation as follows. By assuming the rapid variation of con- 
centrations of clusters larger than the critical size i = 2, their concentrations are approxi- 
mated by those values determined by the given concentrations of the critical cluster as; r$ = 
(k 2 /k 3 )r 2 , r 4 = (k 2 /ki)r 2 + (/x//c 4 a)r|, r 5 = (k 2 /k 5 )r 2 + (/i/k 5 a)r^, r 6 = (k 2 /v)r 2 a+ (jijv)r\, 
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and the corresponding relations for s's, as long as o / 0. In particular, we obtain the relation 
z/rg = A^^Q + A^l- By inserting the relation into the rate equations of the dimers, one obtains 
the reduced rate equations for the dimers as 

r2 = ha 2 - Xr2 + 2/c2^2a + ^r 2 , 

S2 = k\a 2 — \S2 + 2^2520 + ^s 2 .. (10) 

These equations contain essential terms such as the decomposition term — Ar2 (or —XS2) and a 
nonlinear autocatalytic term \ir\ ( or //s^)- In this sense we expect that eq. (10) shows similar 
characteristics to that of eq. (1). For example, if the dissociation is absent, A = 0, then f<i and 
S2 consist only of positive terms and the stationary solution r2 = S2 = is possible only for 
a = T2 = S2 = 0. From the stationary condition for hexamers we obtain r$ = sq = 0, but other 
concentrations rj and Sj with 3 < i < 5 are arbitrary as long as the conservation relation is 
satisfied. No chirality selection is possible, since there are no fixed points, but fixed lines or 
planes. 

We now study the case with dissociation (A > 0). One notices that the fixed points for the 
full rate equation (9) are identical with those obtained by the reduced rate equation (10) by 
way of the construction of the latter. There are racemic fixed points (r 2 = S2) which are not 
of our interest. As for the asymmetric fixed points, it can be easily shown that they should 
satisfy the relations 

A — 2/c2a k\a 2 

T2 + S2 = > 0, r 2 S2 = , 

H n 
c = a + (A + Ba)k 2 (r2 + s 2 ) + (C/a + B)n(r% + sj), (11) 

with A = X^i=2 ^j -1 , = 6/1/, C = X^=4^ rl - The first two relations of eq. (11) yield two 
asymmetric fixed points Si:(r 2 ,+, s 2 ,-) and S2'-{r2-, S2,+) with 



\ — 2k2d I f X — 2k2d\ 2 k\a? 



which is expressed by using a fixed point value of a yet to be determined. For the existence 
of these asymmetric fixed points, a should be in the range X/2(k 2 + \J Jlki) > a > 0. Finally 
the asymmetric fixed points are determined by solving an algebraic equation of o, obtained 
by inserting the first two relations in eq. (11) into the last one. 

Since the conditions for the existence of the asymmetric fixed points are very involved, we 
restrict the discussion here to the case of cq = A;2 = ^3 = &4 = &5 = 1- Then it can be shown 
that the asymmetric fixed points exist for A's in the range 

< A < A u , (13) 

where the upper limit X u is given by the larger root of the quadratic equation 

6/3A 2 + v{^i + 28/3 + 18/3 2 )A - 2/*i/(l + (3) = (14) 
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with (3 = \J [ik\ . In particular, X u ~ 2 for ki ~ and X u ~ /ii//9/3 — > for A;i — > oo . 

We then perform numerical integration of the rate equation (9), in order to observe the 
chiral symmetry breaking in the crystallization. Since the formation of a critical nucleus (= a 
dimer) is expected to be slow, k\ is set as small as k\ = 0.01 compared with other parameters 
as Co = 1, ^2 = ^3 = ^4 = ^5 = 1) A = 0.3, fj, = v = 1. The decomposition rate A = 0.3 is set 
smaller than the upper limit X u = 0.513 calculated from eq. (14), and one of the asymmetric 
fixed points is expected to be realized. This expected behavior is in fact observed in the flow 
diagram shown in Fig. 1 in the phase space of the total number densities of R and S crystals, 
r E = X^=2 i r i an d s s = Ylt=2^ s ii respectively. By starting close to the symmetric mixture 
r E ~ s E , both enantiomeric crystals increase in number, approaching the racemic fixed point, 
but this racemic fixed point is unstable, and eventually one enantiomer overwhelms the other. 
The flow is attracted to one of the chiral fixed points. If one decreases ki, the chiral asymmetry 
represented by the crystalline enantiomeric excess cee = |(r E — s s )/(r s + s E )| increases up to 
almost unity. 
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Fig. 1. Flow diagram in total r's and s's phase space. Symbols represents symmetric and asymmetric 
fixed points. The parameters arc cq = 1, k\ = 0.01, hi = k$ = = k$ = 1, A = 0.3, fi = v = 1. 

To summarize, we proposed a clustering model for the growth of a chiral crystal composed 
of achiral elements, and studied the conditions under which chiral symmetry breaking takes 
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place. The dimers are assumed to be the smallest size possessing chiralities, and at the same 
time of the critical size. When the cluster size increases indefinitely by successive incorporation 
of achiral elements, a racemic mixture results. 

In experiments, chirality selection is achieved by continuous stirring such that large crys- 
tals are fragmented into small clusters and the scattered small clusters induce secondary 
nucleation. In order to simulate the effect of the fragmentation caused by stirring, we intro- 
duce a decomposition process of large clusters into dimers in our crystallization model. If we 
restrict the cluster growth only up to the tetramers, the chiral symmetry is never broken. The 
result is interpreted in the quasi-steady approximation where clusters larger than the critical 
size are assumed to evolve sufficiently rapidly to attain the steady-state concentration. The 
approximation leads to the reduced rate equation of dimers with a linear autocatalysis, which 
is shown to be insufficient for breaking the chiral symmetry, as reported in our previous study. 3 

To equip the the model with the ability to break chiral symmetry dynamically, one needs 
a nonlinear autocatalytic process in the reduced rate equation of the dimers. This is acom- 
plished by incorporating dimers into the other clusters to facilitate cluster growth. The process 
leads effectively to the nonlinear enhancement of the dimer production rate, or a nonlinear 
autocatalytic effect. We know from our previous study that the chiral symmetry can break 
in appropriate situations. 3-5 Numerical integration of the full rate equations confirms this 
expectation. In the present model we included the coagulation of two dimers into a tetramer. 
This process consumes the minimum chiral units or dimers in a nonlinear way, but the frag- 
mentation of a larger cluster, hexamer in our present case, replenishes dimers again into the 
crystallization process. If we incorporate dimer coagulation with clusters of larger sizes as 
R,2+Rj — >Rj+2 and the breaking processes Rj+2 ^R2+Rj, the parameter region for the chi- 
ral symmetry breaking is affected but the qualitative feature remains valid. In the present 
model, the dissociation which corresponds to the back reaction is again found essential, since 
without it the system relaxes to the fixed line in r E — s s phase space, and no precise chirality 
selection takes place. 

Thus far we have assumed that the smallest clusters, dimers, which show chirality are 
at the same time critical clusters which can dissociate, for mathematical simplicity. Actually 
there should be a margin between the two sizes, and the fragmentation of large clusters into 
those with sizes smaller than the critical ones might be the minimal requirement to sustain 
symmetry breaking. This aspect will be studied in the future. 
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